Quenched lattice calculation of semileptonic heavy-light meson form factors 



G.M. de Divitiis"•^ R. Petronzio"•^ N. Tantalo'''^ 

" Umversita di Roma "Tor Vergata", T00133 Rome, Italy 
' INFN sezione di Roma "Tor Vergata", 1-00133 Rome, Raly 
Centro Enrico Fermi, 1-00184 Rome, Italy 

We calculate, in the continuum limit of quenched lattice QCD, the matrix elements of the heavy- 
heavy vector current between heavy-light pseudoscalar meson states. We present the form factors 
for different values of the initial and final meson masses at finite momentum transfer. In particular, 
we calculate the non-perturbative correction to the differential decay rate of the process B Dlu 
including the case of a non-vanishing lepton mass. 



I. INTRODUCTION 



' Semileptonic decays of heavy-light mesons play a central role in the study of flavour physics both on the experimental 
and theoretical sides. The extraction of the Cabibbo-Kobayashi-Maskawa [U, matrix element Vcb, for example, 
: requires the experimental measurement of the decay rate of the process B — > D^*>ii'i and the theoretical calculation 
' of the hadron matrix elements of the flavour changing weak currents. A non-perturbative estimate of the matrix 
^ , elements can be obtained by lattice QCD. Furthermore, within the heavy quark effective theory (HQET) it has been 
O ■ shown 3 that the semileptonic transitions between heavy-light mesons can be parametrized, at leading order of 
, the expansion in the inverse heavy quark mass, in terms of a universal form factor known as Isgur-Wise function. 
■ The Isgur-Wise function is universal in the sense that it describes any semileptonic decay mediated by heavy-heavy 
] weak currents regardless of the flavour of the initial and final heavy quarks and of the spins of the mesons. From 
the phenomenological point of view it is relevant to know the size of the corrections to the Isgur-Wise limit and to 
I establish at which order the heavy quark expansion has to be truncated to produce useful results down to the charm 
'~J^ , mass. 

Matrix elements of the vector heavy-heavy currents between pseudoscalar meson states are parametrized in terms of 
(— I ■ two form factors. In the case of the light leptons £ = e, ii, the differential decay rate of the process B Dive is 
I— I, proportional to the square of a particular linear combination of the two, G^^^. The BaBar and Belle collaborations 
have already measured [1, 0] the branching ratios of the processes Ti/r and a future measurement of the 

differential decay rate will make possible to extract Vcb also from this channel. In this case a separate knowledge of 
the form factors is required, both in the Standard Model and in its minimal extensions (see for example refs. [6|,[7[). In 
ref. [sj we have already shown our final results for G^~*^ by focusing on their phenomenological implications without 
giving all the details of the calculation. Here we present detailed results of continuum and chiral extrapolations and 
' separate estimates of the two independent form factors for several values of the initial and final heavy quark masses 
J — , together with an analysis of the infinite heavy quark mass limit. In addition, we make a prediction for the ratio of 
the differential decay rates of the processes B — > Dlvg with i — t and ^ = e, 

The simulation of relativistic heavy quarks with masses ranging from the physical b mass down to the physical c mass 
has been performed by using the step scaling method (SSM) [9||, already applied successfully to the determination 
of heavy quark masses and heavy- light meson decay constants [10, [HI, [iJ]- The SSM allows to reconcile large quark 
masses with adequate lattice resolution and large physical volumes. The two form factors have been calculated for 
different values of the momentum transfer by making use of flavour twisted boundary conditions [l3] , that shift the 
d ^ discretized set of lattice momenta by an arbitrary amount (see also fisl [l6{). 

The plan of the paper is as follows. In section [TT] we introduce the form factors in the continuum theory and re- 
derive the Luke's theorem . In sections HII] and |TV] we set up the lattice notation and describe the calculation. 
In sections |V] we discuss the results at finite volumes while in section IVII we show our final results. We draw our 
conclusions in section fVIll 
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II. FORM FACTORS 



The semileptonic decay of a pseudoscalar meson into another pseudoscalar meson is mediated by the vector part of 
the weak V — A current. The corresponding matrix element can be parametrized in terms of two form factors, 

{Mf\ \M.) = {p.+Pfr /r^ + ip^-Pff f--^^ 

or, equivalently, 

where Vij — pij / Mij are the 4- velocities of the mesons. The relations between the h^^'^ and the parametriza- 
tions are given by 

^.^/ _ (M, + M^)/r^ + (M,-M^)/- ^ 
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2y/M~Mf 



^ 2^M,Mf 

In the rest of this paper we work in the h^^^ parametrization that, as it will emerge from the discussion below, is 
more convenient for the study of the dependence of the form factors upon the masses of the initial and final heavy 
quarks. 

The form factors depend upon the masses of the parent and daughter particles and upon w = Vf ■ Vi 

h'^^'iw) = h±{w,M,,Mf), 

Time reversal and hermiticity imply that h^_^'^ and h^^^ are real. Furthermore they imply that h^^'^ is even under 
the interchange of the initial and final states while h^Z*^ is odd, 

h+{w,Nh,Mf) = h+{w,Mf,M{), h.{w,M„Mf) = - h^{w,Mf,M,) (2) 

In eq. ([T]) one can consider the limit in which both meson masses go to infinity at fixed 4- velocity; the left hand side is 
well defined in this limit and, consequently, also the form factors. It is thus legitimate to make a change of variables 
from the meson masses to the parameters e+ and e_, defined as 

11 _ 1 1 

Expressed as functions of the new variables, h^^^ = h±{w, e+, £_) are well defined at = and e_ = and can be 
expanded in power series around these points. The symmetry properties of eq. ([2]) force the odd(cven) powers of £_ 
to vanish into the expansion 

h^{w,s^,s^) = /HKO,0)+e+ +- + - ^_ + ... 

/.-Ks+,£_) = e + ^^^^ + ... (3) 

In the elastic case, when the initial and final mesons coincide, h^_Z*^ vanishes and the vector current is conserved. The 
conservation of the vector current implies that h^_^^{w = 1) = 1. This condition, inserted in the previous equations, 
translates into a condition on the derivatives of /i!^* with respect to e+ 

9"/t+(w= 1,0,0) _ 

del 



We thus expect that, for values of w ~ 1 the corrections proportional to e+ will be rather small while the ones 
proportional to e_, not constrained by the vector symmetry, can play a role also at zero recoil. This expectation is 
confirmed by our numerical results (see section fyi]) . 

The discussion above is a re-derivation of the "Luke's theorem" [17]. The theorem, originally derived by using HQET 
arguments, states that h^_^'^ it is not affected by first order corrections at zero recoil. In our language 

I, / . ^ . el d^h+{w^ 1,0,0) 

h+{w = l,e+,e-) = 1 + — + (4) 

Eqs.^]) confirm the analysis of the subleading corrections to the form factors that has been carried out by the authors of 
ref. [18] within HQET^. The additional symmetries of the static theory imply relations among the coefficients appearing 
in eqs. ^ and the corresponding ones arising in the case of vector-pseudoscalar and vector-vector transitions. In 
particular, h^_^^ is proportional at leading order to the Isgur-Wise function 0]. 

In ref. Q we have shown the results concerning the form factor G^^^{w) that enters into the semileptonic decay 
rate of a i? meson into a D meson in the approximation of massless leptons £ — e, ^, 

Kw < — ^ ^ 

~ ~ 2MbMd 

This form factor is related to h^^^ {w) and h^T^ {w) by 

In the case I = t the mass of the lepton cannot be neglected and the differential decay rate is given by [1, [l^ 



dw dw \ t{w) J \\ 2t{w) J 2t{w) w 



rr^-ri~^ r=— -, t{w) = l + r -2rw, 
Mb Mb 



Ml + Ml - ml 



2MbMd 

In this work we provide an estimate of the function A^^^(u') appearing in the previous relations, including values 
at It; > 1. Its expression in terms of hJ'T*^ {w) and h^^'^ (w) is given by 
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J i^ff . w-1 ' 

/i I (w) — n (w) 

1+r + ^ ' w+1 " ^ ' 
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1-r w-1 h'Z"^{w) \ ( 1 - r h'T^ {w) 



^ Some care is needed when eqs. l(3)l are compared with the corresponding results of ref. 18|. Indeed eqs. ((S} are the result of a Taylor 
expansion and the coefficients do not depend upon and, consequently, upon the meson masses. Eqs. (Bl) of ref. [Tsll are expansions 

in inverse powers of the quark masses that depend upon the renormalization scale as well as the coefficients. This dependence cancels 
at any given order. 



In the elastic case A*^'^(w) vanishes identically and, in the approximation in which h^Z*^ {w) is much smaller than 
K^'^ (w), it is very well approximated by its static limit 

A^H.i-I, r = ^ (6) 

l + r Mi 



III. LATTICE OBSERVABLES 



We have carried out the calculation within the 0{a) improved Schrodinger Functional formalism [20|, |21| with T ~ 
2L and vanishing background fields. Physical units have been set by using the Sommer's scale and fixing ro — 
0.5 fm [HjlHjHj- In order to set the notations, we introduce the following source operators 



where s and r are flavour indexes while and C are boundary fields at xq — and xq — T respectively. The bulk 
operators are defined according to 

= ^,(x)757Vr(x), Psr{x) = M^h^M^^) -^lix) = ^"r(^) + aCA^^^Psr{x) 

V^Ax) = Mx)rMx), = ^s{x)rYMx) V]^,{x) = Vt^x) + acv^^^^Tr(x) 



The improvement coeflicient ca has been computed non-perturbatively in ref . [25| . Regarding cv , we have used the 
perturbative result from ref. |26| but its actual value influences our results at the level of a few per mille. 

The quark masses have been defined through the PCAC relation. We have calculated the following correlation 
functions 



f^Ax,)=-Y.^OrsAl{x)) 



ffAxo)=^Y.(OrsPsrix)) 



and defined 
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/4. + acAdodgf^r 
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where a is the lattice spacing, kr is the hopping parameter of the r quark and fee is the critical value of the hopping 
parameter. The renormalization group invariant (RGI) quark masses have been obtained by the following relation 



= Zm [1 + [hA - bp) amX"^'] 



AW I 



(7) 



The combination hA — bp of the improvement coefficients of the axial current and pseudoscalar density has been 
computed non-perturbatively in [27l . [29l | . The factor Zm is known with very high precision in a range of inverse bare 
couplings that does not cover all the values of (3 used in our simulations. We have used the results reported in table 6 
of ref. 28] to parametrize Zm in the enlarged range of j3 values [5.9, 7.6]. 

In order to define on the lattice the matrix elements of the vector current between pseudoscalar meson states, we need 
to introduce other two correlation functions. 



/i(xo,Pf) = -5](0,M°;(x)) 



where i and / are the heavy fiavour indexes and / is the light one. The external momenta have been set by using 
fiavour twisted b.c. for the heavy fiavours; in particular we have used 



i^ijix + IL) 



'i^^,fix), 



Pi 



L 



ki G 



and ordinary periodic b.c. in the other spatial directions and for the hght quarks. We have worked in the Lorentz 
frame in which the parent particle is at rest (pi = 0). In this frame w is simply expressed in terms of the ratio between 
the energy and the mass of the daughter particle w = Ef/Mf. The matrix elements of have been defined by the 
following ratios 

{Vn7/ ^ {Mf\ \M..)ni ^ 2VM^^ , ■Ff_^(r/2; 0, pf) 

/.F0^.(T/2;0,0).Fj?^^(T/2;pf,pf) 

that become the physical matrix elements in large volumes where single state dominance is a good approximation. 
An alternative definition of the matrix elements {D2), which reduces to the previous one {Dl) in the infinite volume 
and at zero lattice spacing, can be obtained by considering 



_^ VM-Eff^iT/2,0) .ff_/r/2;0,pf) 

/M7/_^(T/2,pf) Jj^o_^^(T/2; 0, 0).F0^/r/2; 0, 0) 



In eqs. ([HI) and Q the renormalization factors Zy and Za cancel in the ratios together with the factors containing 
the improvement coefficients hy and hA- 

By calculating the following ratio 



^/%(r/2;0,pf) {Mf\VO\Mf) 



w - 



we have defined w, as well as meson masses and energies, entirely in terms of three point correlation functions. This 
definition of w is noisier than the one that can be obtained in terms of ratios of two point correlation functions; 
however it leads to exact vector current conservation when Mf = Mi and reduces the final statistical error on the 
form factors. The two definitions of the matrix elements lead to two definitions of the form factors that, in terms of 
{V")d and {V^)d, are expressed by 



The last two equations are not defined at w = 1; this is due to the second term in the parenthesis of eq. (fTTj) and p2|) 
that we extrapolate at zero recoil before calculating h^Z*-^ [w = 1) and G^^^ [w = 1). 



IV. THE STEP SCALING METHOD 

The SSM has been introduced to cope with two-scale problems in lattice QCD. In the calculation of heavy-light meson 
properties the two scales are the mass of the heavy quarks (fo,c) and the mass of the light quarks {u,d,s). Here we 
consider the generic form factor F^^^f = {h";^ ^ , h^T ^ , G'-^ f } as a function of w, the volume and fix the meson 
states by the corresponding heavy and light RGI quark masses that, being extracted by the lattice version of the 
PC AC relation, arc not affected by finite volume effects (see eq. [7]). 

The first step of the finite volume recursion consists in calculating the observable F"^^^ {w; Lq) on a small volume, 
Lq, which is chosen to accommodate the dynamics of heavy quarks with masses ranging from the physical value of 
the charm mass up to the mass of the bottom. As in our previous work we fixed Lq — 0.4 fm. We have simulated 
five different heavy quark masses niij — {■m\,m\,m'l,mf^,ml^}, five different momenta 6i = {6l,9f,9l,0f,6l} (see 
eq. ([5])) and three light quark masses mi = {mj,mf,mf}. 
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FIG. 1: Extrapolations to the continuum limit of G^^^ (w). The data correspond to m; — nis, to the definition Dl and to 
the data sets LqA, LqB and LqC. 



A first effect of finite volume is taken into account by evolving the results from Lq to Li = 0.8 fm through the factor 



F^^f{w;Lo) 

computed for each value of w and for each value of the light quark mass. The crucial point is that the step scaling 
functions are calculated by simulating heavy quark masses smaller than the &-quark mass. The step scaling functions 
at rrii ~ rUb and ruf rric are obtained by directly simulating uif both on Lq and on Li and by a smooth extrapolation 
in l/rrii. 

Extrapolating the step scaling functions is more advantageous than extrapolating the form factors. This can be easily 
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FIG. 2; Light quark mass dependence of h^~^^ {w\ Lq) (left) and of h^^^ {w; Lo) (right). The different sets of points correspond 
to different values of rrii ranging from about tUs to about tUs/A. The data are in the continuum limit and correspond to the 
definition Dl. 



understood by relying on HQET expectations (see also eq. ^ 
a'^f{w;Lo,Li) 
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In the previous relations the superscripts in parenthesis, (n), mark the order of the expansion in the inverse heavy 
quark mass. The subleading correction to the step scaling functions is the difference of two terms and vanishes in the 
infinite volume, a^^^^^ {w; Lq, Li) — F^^^^^^ {w; Li) — F'^^')^^ {w; Lq), becoming smaller and smaller as the volume is 
increased. This matches the general idea that finite volume effects, measured by the cr's, are almost insensitive to the 
high energy scale. 

We also compute the step scaling functions of the elastic form factors /I'p* at rrii — by extrapolating the cor- 
responding results from smaller heavy quark masses. Also in this case the cr's are expected to be almost flat with 
respect to I/to^. 

In order to remove the residual finite volume effects we iterate the procedure described above once more passing from 
ii to L2 = 1-2 fm. Our final results arc obtained from 



F'^f{w-L2) 



F'-f{w;Lo) a'-fiw;Lo,Li) a'^^ {w; L^, L2) 



(13) 



V. FINITE VOLUME RESULTS 



A. small volume 



The small volume Lq — 0.4 fm has been simulated by using three different values of the lattice spacing (see table HB. 
The small physical extent of the volume allowed us to simulate relativistic heavy quarks with masses ranging from 
around down to rric- We have computed the form factors /i^^, /i!r*'^ and G^~*^ for all the combinations of heavy 
and light quark masses and for five different values of the momentum transfer. 

In figured] we show the continuum extrapolations of G^^^ (w). The points in this figure correspond to to/ = but 
similar figures can be obtained for the other values of the light and heavy quark masses and for the other form factors. 
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FIG. 3: Continuum extrapolation of a]^''{'w = 1.1; Lq, Li) (left) and (T!jr*'(w = 1.1; Li, L2) (right) at the heaviest values of the 
heavy quark masses (m; ~ m(,/4 and rrii ~ mi,/2 respectively). The data correspond to mi — nis, to the definition Dl and to 
the data sets LiA/Loa, LiB/Lob (left) and L2A/L\a, L2B / Lib (right). 
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FIG. 4: Step scaling functions of h'J^'^ (left) and /i!r*'^ (right) as functions of l/rui for the first evolution step (from Lq to Li). 
The black vertical lines represent the physical points rui — rric and rui — mi,. The data are in the continuum and chiral limits 
and correspond to the definition Dl. 



In figure [2] we show h^~^^ (left) and h^~^^ (right) as functions of w for the three different values of light quark 
masses that we have simulated (ranging from about rus to ms/4). As we have anticipated in ref. @, we find that 
the -F's behave as constants with respect to mi within the statistical errors. This happens for each combination of 
heavy quark masses and for each value of the lattice spacing. Nevertheless we make a linear extrapolation to reach 
the chiral limit; the resulting error largely accounts for the systematics due to these extrapolations. In the following 
our results include the mild chiral extrapolation. 



B. steps 



The parameters of the simulations of the evolution steps are given in tables IIVI and IVII . We have been simulating at 
two different lattice spacings by limiting the maximum value of the heavy quark mass to rrii ~ TOf,/2 for the first step 
and to rrii — mi,/ 4: for the second. In figure [3] we show the dependence upon the lattice spacing of the step scaling 
functions a'^_^-^{w = l.l;Lo,Li) (left) and cr!p'^(w = l.l;Li,L2) (right) in the worst case (largest value of heavy 
quark masses). In general our results are consistent with a scaling regime within a few per mille accuracy and the 
continuum step scaling functions of tables |V] and IVIII have been obtained by averaging the results at the two lattice 
spacings. 

In figure m we can test our hypothesis on the low sensitivity of the step scaling functions upon the high energy scale. 
The figure shows the step scaling functions of the form factors h^^'^ (left) and /i!!*"^ (right) as functions of 1/mi. In 
both cases the dependence upon rrii is hardly appreciable and in the case of /i*^"^ the cr's are very close to one while 
affected by stronger finite volume effects. We obtain the values at mi — mi, by linear fits. 
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FIG. 5; Step scaling functions of h^_^^ (left) and h^S*^ (right) at fixed m/ as functions of l/rrii for the second evolution step 
(from Li to 1/2)- The black vertical lines represent the physical points rrii — rric and rrii — mt- The data are in the continuum 
and chiral limits and correspond to the definition Dl. 



1.2 
1.1 
1 

0.9 
0.8 
0.7 



* S 2 




9 










e 


'di a 






D2 - » : 







1.2 
1.1 

1 

0.9 
0.8 
0.7 













i 


'di O 




D2 : » : 





1.1 1.15 1.2 1.25 



1.15 1.2 



FIG. 6: Comparison of the two definitions of h^^^ {w; L) at Lo — 0.4 fm (left) and at L2 — 1.2 fm (right). The data are in 
the continuum and chiral limits. 

In figure Owe plot the same quantities as in figure S] for the second evolution step (from Li to L2, see table IVl| ) . Also 
in this case the step scaling functions depend very smoothly upon I /mi. 



C. consistency checks 



In this section we illustrate the results of two checks that we have done in order to convince ourselves on the consistency 
of the step scaling procedure. As already discussed in sec. IIIH we have used two different definitions of the matrix 
elements and, consequently, of each form factor. In figure [S] we show the comparison of h^^^ [w; L) at Lq — 0.4 fm 
(left) and at L2 — 1.2 fm (right). We see that the results, while differing at finite volume, converge to common values 
after the step scaling procedure. This makes us confident of a correct accounting of finite volume effects. 

A second check of the whole procedure, and in particular of the continuum limit of the step scaling functions, can 
be obtained by considering the elastic form factor h^^^ {w; L) at fixed w as a function of L. The point is that the 
charm quark mass has been simulated directly on each physical volume and, in particular, on the biggest one. In 
figure [7] we fix w = 1.05 (left) and w = 1.10 (right) and see that the step scaling recursion (black points) converge 
to the result obtained directly at L2 = 1.2 fm (red points, slightly displaced to help the eye) making us confident of 
a correct accounting of the cutoff effects and, in particular, of a correct estimate of the error on the continuum step 
scaling functions. 



Our final results are obtained by averaging over the two definitions and by combining in quadrature statistical errors 
with the systematic ones that we estimate from the dispersion between Dl and D2. 
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FIG. 7; = 1.05;!/) (left) and h'^^^iw — 1.10;!/) (right) as functions of the volume. The black points have been 

obtained through the step scaling recursion while the red points (slightly displaced on the a;-axis to help the eye) are the result 
of a direct simulation on the biggest volume. The data are in the continuum and chiral limits and correspond to the definition 
Dl. 
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FIG. 8: In the left plot it is shown h^^^ {w) in the infinite volume limit as a function of w for different values of the final heavy 



quark mass. The right plot shows (w) for the same combinations of heavy quark masses. 



VI. FINAL RESULTS 



In this section we discuss our final results in the continuum, chiral and infinite volume limits (table |T]) . In order to 
establish the onset of the static limit approximation we plot in figure [S] the form factor h^~^^ {w) as a function of w 
for different values of the final heavy quark mass. The right plot shows h^~^^ {w) for the same combinations of heavy 
quark masses. We see that the corrections to the static limits of both h^^^ {w) and h^^^ (w) are of the order of 2% 
at the charm mass. For heavy quark masses bigger than mf,/2 the corrections are almost negligible (below 1%). 

Eqs. ^ and ^ predict that the convergence toward the static limit is faster in the case of the elastic form factors with 
respect to the ones having mi > ruf. This happens because near the point at zero recoil the subleading corrections to 
li^^ (w) are proportional to the square of the difference of the initial and final meson masses. Figure [HI clearly shows 
that this happens in practice. Indeed, in the left plot we see that the elastic form factor h^^^{w) is much closer 
to the static limit (very well approximated by h^^^{w)) with respect to the form factor h^^-^{w), the one relevant 
into the calculation of Vcb shown in figure \8\ In the right plot of figure [9] we show how well eq. (jl]) is approximated 
by our numerical data. The fit is performed on the slope while the intercept is fixed to one. 

The QCD form factor KZ^Hw) is related to the renormalization group invariant HQET Isgur-Wise function, f (w), by 
the following relation [30, l3l| 

h''^^{w)= 1 + /3+(mj, m/; uj) + 7+(mj, m/; w) + O(TOj^y) ^(w) 
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and (w): in the range 1 < w < 1.05 the two elastic form factors are indistinguishable 
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within the quoted errors while h^~*^{w) in fig. |5] shows appreciable corrections from the Isgur-Wise limit, in particular at zero 
recoil. The right plot shows /i^^ at zero recoil {w = 1) as a function of e?. (actually {l/nii — l/nif)'^ oc el, rriij being the 
RGI heavy quark masses): the solid line has been obtained by fitting the data according to eq. 
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FIG. 10: The figure shows the function A^^^(w) in the chiral, continuum, and infinite volume limits. The solid line correspond 
to the static limit result, {Mb — Mo)/ {Mb + Md), and has been drawn by using the experimental determinations of the meson 
masses. 



where the 7+ term accounts for non-perturbative power corrections proportional to the inverse of the quark masses 
while the /?+ term accounts for perturbative radiative corrections. In the case of the elastic form factor h[^^{w) at 
the highest value of the simulated heavy quark masses, i.e. the bottom quark mass, power corrections are completely 
negligible in our data as clearly emerges from figures [S] and [H 

7+(TOh,TOfc; w) = 



/if^^H = [l + /3+(mfc,mfa;u.)]^H 

The function P^{mb,mb',w) depends logarithmically upon the bottom mass through as{mii) and vanishes at zero 
recoil where the renormalized Isgur-Wise function is identically equal to one like the relativistic QCD elastic form 
factor. These terms are of the percent order and their logarithmic dependence upon mf, cannot be extrapolated away 
from our data. Nevertheless, in order to get the HQET Isgur-Wise function our non-perturbative results for h^~^^{w) 
(given in table H} can be further corrected by hand through the perturbative l3+{mb,rnb]w) given in ref. [sO] at next 



i f 


w 


G 




h- 


A 




1.00 


1 r\r\r\ / r\r\\ 

1.000(00) 


"1 r\r\r\ ( r\r\\ 

1.000(00) 








1.03 


0.971(07) 


0.971(07) 






D 


1.05 
1.10 
1.20 


0.955(06) 
0.916(09) 
0.828(20) 


0.955(06) 
0.916(09) 
0.828(20) 








1.00 


1 r\r\r\ I r\r\\ 

1.000(00) 


"1 r\r\r\((\r\\ 

1.000(00) 








1.03 


0.974(07) 


0.974(07) 






B ^ B 


1.05 
1.10 
1.20 


0.952(07) 
0.903(16) 
0.794(34) 


0.952(07) 
0.903(16) 
0.794(34) 








1.00 


1.026(17) 


1.017(03) 


-0.011(23) 


0.466(26) 




1.03 


1 r\r\ 1 / 1 r\\ 

1.001(19) 


0.986(08) 


-0.018(19) 


0.465(25) 


B ^ D 


1.05 


0.987(15) 


0.970(07) 


-0.023(16) 


f\ A f A / C\ A\ 

0.464(24) 




1.10 


0.943(11) 


0.928(10) 


-0.024(12) 


0.463(24) 




1.20 


0.853(21) 


0.835(21) 


-0.018(13) 


0.463(23) 



TABLE I: Physical results. Average of the two definitions Dl and D2. 



to leading order^. 

Finally, we show in figure fTUl our best result for the function A^^^(w) that enters in the decay rate of the process 
B DtVt (see discussion at the end of section ITU) . ^^^^{w) does not show any significant dependence upon w and 
is very well approximated by its static limit (see eq. [S]). These findings represent a prediction that can be confirmed 
by a future measurement of the differential decay rate of the process B Drvr- Indeed, the function A-^^^(u') 
can be extracted experimentally by the ratio dV^'^^'^'^^ / dT^'^^'^'^'^^^'^'' i^ that does not depend upon the CKM matrix 
element. On the other hand, the knowledge of A^~*^(u') is required in order to perform lepton-fiavour universality 
checks on the extraction of Vcb- 



VII. CONCLUSIONS 

We have performed the calculation of the form factors that parametrize semileptonic transitions among pseudoscalar 
heavy-light mesons and made a prediction for the ratio dV^^^'^"^ / dT^^^'''^''^^'^'' i' . In view of a future measurement 
of the differential decay rate of the process B — > DrVr, our results will allow to perform lepton-flavour universality 
checks on the extraction of Vcb- 

The form factors have been obtained with a relative accuracy of the order of a few percent allowing to establish the 
range of validity of the heavy quark effective theory for these quantities. In particular we have obtained a check of 
the predictions of the Luke's theorem that we re-derived. The corrections to the static limit are very small already in 
the case of the elastic form factor h^^^ and negligible in the case of h^^^ . We have also established the accuracy 
of the static approximation to the form factors of the decay B Div which is of the order of 2-3% at zero recoil and 
reaches about 7% at ui = 1.2 where becomes definitely inadequate for precise phenomenological applications. 

Our results have been obtained within the quenched approximation and further calculations will be needed to asses 
the corrections due to unquenching. On the other hand, the accuracy reached in the quenched case demonstrates the 
feasibility and the opportunity of repeating the present calculation in the unquenched theory. Indeed, the recursive 
matching process can be extended to the sea quark masses that, alternatively, can be kept to their physical values 
if the Schrodinger Functional formalism is used. Moreover, flavour twisted boundary conditions can be used for 
heavy valence quarks also in the Nf = 3 unquenched theory. The real case will further differ by the heavy flavour 
determinants that can be accounted for by a perturbative expansion in the hopping parameter. 



^ for a recent lattice calculation of the Isgur-Wise function see ref. [32l | 
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i TxL^ iVc„/g k 

LqA 7.300 48 X 24^ 277 0.124176 0.000000 

0.124844 0.953456 

0.128440 1.201080 

0.129224 2.042983 

0.131950 2.573569 
0.134041 
0.134098 

0.134155 

LqB 7.151 40 X 20^ 224 0.122666 0.000000 

0.123437 0.953456 

0.127605 1.201079 

0.131511 1.719170 

0.131686 2.488490 
0.134277 
0.134350 

0.134422 

LoC 6.963 32 x 16^ 403 0.120081 0.000000 

0.120988 0.953456 

0.126050 1.201079 

0.131082 1.719170 

0.131314 2.488490 
0.134526 
0.134614 
0.134702 



TABLE II: Table of lattice simulations of the small volume. 





w 


G 


h+ 


h- 




1.00 


1.000(00) 


1.000(00) 






1.03 


0.979(02) 


0.979(02) 




D-^ D 


1.05 


0.967(03) 


0.967(03) 






1.10 


0.942(04) 


0.942(04) 






1.20 


0.894(06) 


0.894(06) 






1.00 


1.000(00) 


1.000(00) 






1.03 


0.980(02) 


0.980(02) 






1.05 


0.969(03) 


0.969(03) 






1.10 


0.940(05) 


0.940(05) 






1.20 


0.882(09) 


0.882(09) 






1.00 


1.025(17) 


1.013(03) 


-0.020(37) 




1.03 


1.009(14) 


0.992(03) 


-0.032(29) 


B 


1.05 


1.000(13) 


0.980(04) 


-0.040(23) 




1.10 


0.976(11) 


0.953(04) 


-0.048(17) 




1.20 


0.929(09) 


0.903(06) 


-0.053(12) 



TABLE III: Small volume results, Lo = 0.4 ftn. Results corresponding to the definition Dl. 





B 


T X 


Nrr,fa 
^ cnjg 


k 


e 


L/Qd 


6.737 


24 X 12^ 


608 


0.12490 


0.000000 










0.12600 


0.953456 










0.12770 


1.201080 










0.12979 


1.719172 










0.13015 


2.488491 










0.13430 












0.13460 












0.13490 




Lob 


6.420 


16 X 8^ 


800 


0.120674 


0.000000 










0.122220 


0.953456 










0.124410 


1.201079 










0.127985 


1.719172 










0.128066 


2.488491 










0.134304 












0.134770 












0.135221 






6.737 


48 X 24^ 


260 


0.12490 


0.000000 










0.12600 


2.042983 










0.12770 


2.573569 










0.12979 


3.438340 










0.13015 


4.976980 










0.13430 












0.13460 












0.13490 




LiB 


6.420 


32 X 16^ 


350 


0.120674 


0.000000 










0.122220 


2.042983 










0.124410 


2.573569 










0.127985 


3.438340 










0.128066 


4.976980 










0.134304 












0.134770 












0.135221 





TABLE IV: Table of lattice simulations of the first step. 





w 


era 




a- 




1.00 


1 r\r\r\ ( r\r\\ 

1.000(00) 


■1 r\r\c\ ( r\tW 

1.000(00) 






1.03 


0.999(01) 


0.999(01) 




D D 


1.05 


0.997(02) 


0.997(02) 






1.10 


0.993(03) 


0.993(03) 






1.20 


0.985(05) 


0.985(05) 






1.00 


■1 r\r\r\ / c\c\\ 

1.000(00) 


T r\r\r\ f r\r\\ 

1.000(00) 






1.03 


0.997(02) 


0.997(02) 




B ^ B 


1.05 


0.996(02) 


0.996(02) 






1.10 


0.991(04) 


0.991(04) 






1.20 


0.981(09) 


0.981(09) 






1.00 


■1 c\c\r\ f r\c\\ 

1.002(02) 


1.003(01) 


0.89(16) 




1.03 


0.999(03) 


T r\r\r\ / 

1.000(02) 


0.86(13) 


B —> D 


1.05 


0.996(04) 


0.998(02) 


0.83(11) 




1.10 


0.991(04) 


0.993(03) 


0.77(08) 




1.20 


0.980(05) 


0.983(05) 


0.72(08) 



TABLE V: First step, from Lo = 0.4 fm to L\ = 0.8 fm. Results corresponding to the definition D\. 








T X 


Nrr,fa 
^ cnjg 


k 


e 




6.420 


32 X 16^ 


360 


0.126600 


0.000000 










0.127400 


1.603930 










0.128030 


2.080840 










0.128650 


2.978423 










0.129500 


4.311249 










0.134304 












0.134770 












0.135221 




Lib 


5.960 


16 X 8^ 


480 


0.118128 


0.000000 










0.119112 


1.603930 










0.120112 


2.080840 










0.121012 


2.978423 










0.122513 


4.311249 










0.131457 












0.132335 












0.133226 




L2A 


6.420 


48 X 24^ 


250 


0.126600 


0.000000 










0.127400 


2.405895 










0.128030 


3.121260 










0.128650 


4.467634 










0.129500 


6.200000 










0.134304 












0.134770 












0.135221 




L2B 


5.960 


24 X 12^ 


592 


0.118128 


0.000000 










0.119112 


2.405895 










0.120112 


3.121260 










0.121012 


4.467634 










0.122513 


6.200000 










0.131457 












0.132335 












0.133226 





TABLE VI: Table of lattice simulations of the second step. 





w 


era 


0^+ 






1.00 


1 r\r\r\ f r\r\\ 

1.000(00) 


■1 r\r\c\ ( r\tW 

1.000(00) 






1.03 


U.99j(01) 


0.993(01) 




D D 


1.05 


0.988(02) 


0.988(02) 






1.10 


0.973(05) 


0.973(05) 






1.20 


0.921(16) 


0.921(16) 






1.00 


■1 c\r\r\ { r\r\\ 

1.000(00) 


T AAA/AA\ 

1.000(00) 






1.03 


0.992(02) 


0.992(02) 




B ^ B 


1.05 


0.986(03) 


0.986(03) 






1.10 


0.961(10) 


0.961(10) 






1.20 


0.890(24) 


0.890(24) 






1.00 


"1 r\c\r\ ( c\A \ 

1.000(01) 


"1 AAA/ A1 \ 

1.000(01) 


0.63(32) 




1.03 


0.992(02) 


0.993(02) 


0.65(25) 


B —> D 


1.05 


0.987(02) 


r\ aoo/aa\ 

0.988(02) 


0.67(21) 




1.10 


0.972(05) 


0.972(05) 


0.65(21) 




1.20 


0.921(14) 


0.921(14) 


0.45(33) 



TABLE VII: Second step, from L\ = 0.8 fm to L2 = 1.2 fm. Results corresponding to the definition Dl. 



